Week #13 È بداعي‎ 
Frequency Response الملل‎ 


The Concept of Frequency Response 


T 


- In the steady state, sinusoidal inputs to a linear system generate sinusoidal responses of 
the same frequency. Even though these responses are of the same frequency as the input, they 
differ in amplitude and phase angle from the input. 


- Sinusoids can be represented as complex numbers called phasors. The magnitude of the 
complex number is the amplitude of the sinusoid, and the angle of the complex number is the 
phase angle of the sinusoid. 


Mı cos (wt + $) represented as M,Z¢, 


Consider the mechanical system of figure below. If the input force, f(t), is sinusoidal, the 
Steady-state output response, x(t), of the system is also sinusoidal and at the same frequency as 
the input. 
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2- Polar plot 


MATLAB code 


Syms s w 
G= 1/(s+2); 
pretty (G) 


s= (۷۷ ع‎ 

G= subs (G,s) ; 
w—[0.1:100]; 
G= subs (G,w) ; 


M =abs (G) 
Fi =rad2deg (angle )6( ( 


M= 20*10g10 (M) 
subplot(211),semilogx(w,M) 
xlabel('Frequency (rad/s') 
ylabel('20 Log (M) ') 

title ('Magnitude') 

grid 


subplot (212) , semilogx(w, Fi) 


xlabel('Frequency (rad/s') 
ylabel ) ۲ Phase (degree) ') 
title('Phase') 

grid 

figure 


` compass (G) 
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using log-frequency in rad/s as the ordinate, 


log-magnitude and the phase, 


b- Make plots of the 


G(s) = (s +2)(s + 4) 
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«sions for the magnitude and phase responses of 


Find analytical expre‏ ۔و 


EX 





f'Asymptotic Approximations: Bode Plots 


The log-magnitude and phase frequency response curves as functions of log «arc called Bode 
plots or Bode diagrams. Sketching Bode plots can be simplified because they can be 
approximated as a Sequence of straight lines. Straight-line approximations simplify the 
evaluation of the magnitude-and phase frequency response. 

Bode Plots for G(s) = (s + a) 


Consider a function, G(s) = (s + a), for which we want to sketch separate logarith- 
mic magnitude and phase response plots. Letting s = jw, we have 


Oi ja) = (jw +a) = ره‎ + i) 
At low frequencies when w approaches zero, 
G(jo) = a 
The magnitude response in dB is 
20log M = 20 log a 
where M = |G(je»)| and is a constant. . | 
At high frequencies where w > a. 

G(jo) = (A) = a) 90° = w ”او ے‎ 
The magnitude response in dB is 

20 log M = 20 log 4 + 20 log ^ = 20 loge 
where a > w > oo. Notice from the middle term that the high-frequency approxi- 


mation is equal to the low-frequency approximation when o =a, and increases 
for w> a. 


wee s —‏ وس وی > رک ی ی س ی ج سے د کے سس - 


If we plot dB, 20 log M, against log w, 
20108 M = 20logo 
y = 20x 


wheres y = 201og M, and x = log o» The line has a slope of 20 when plotted as dB vs. log w. 


© Since each doubling of frequency causes 20 log w to increase by 6 dB, the line 
rises at an equivalent slope of 6 dB/octave, where an octave is a doubling of 
frequency. This rise begins at م‎ = 4, where the low-frequency approximation equals 
the high-frequency approximation. IT 
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Many times it is convenient to draw the line over a decade rather than an octave, where a | 
decade is 10 times the initial frequency. Over one decade, 20 log « increases by 20 dB. Thus, 4 
slope of 6 dB/octave is equivalent to a slope of 20 dB/decade. 
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Bode plots of (s + a): 4. magnitude plot; b. phase plot 
Note 


It is often convenient to normalize the magnitude and scale the frequency so that the log-magnitude 
plot will be 0 dB at a break frequency of unity. 

To normalize (s + a), we factor out the quantity a and form a[(s/a) + 1]. The trequency is scaled by 
defining a new frequency variable, s, = s/a. Then the magnitude is divided by the quantity a to yield 

0 dB at the break frequency. Hence, uic normaiized and scaled function is (s1 +1). To obtain the 
original frequency response, the magnitude and frequency are multiplied by the quantity a. 
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Bode magnitude plot: slope contribution from each pole and zero 
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MATLAB code 


Gi-tf(i,poly([ 0 .5 -41)); 
G2-tf([1 5],poly([ وس‎ -41)) ; 
Gort (polv([-3 -5]) ,Poly([ 0 
zpk (G1) -2 -4])); 
zpk (G2) 
zpk (G3) 


‘subplot (131) , 


bode (G1): ; 
title( 'G1z1/ (s ); grid : 
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subplot (133) " bode (G3) ; grid 
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Stability analysis in a Bode diagram 


It is most common to use a Bode diagram for frequency response based stability analysis 
closed loop systems. The Nyquist's Stability Criterion says: The closed loop system وز‎ 
marginally stable if the Nyquist curve (of L) goes through the critical point, which is the Poin, 


(-1, 0). 


But where is the critical point in the Bode diagram? The critical point has phase (angle) - 0۰ 
and amplitude 1 = OdB. The critical point therefore constitutes two lines in a Bode diagram:~ 
The OdB line in the amplitude diagram and the -180 » line in the phase diagram. 


[dB] 






0 dB 






£ 





(0 
(logarithmic) 






[degrees] 


۰ 
eee Lee RN ئا‎ tere, 


- 180 


PM = ZI (jue -)- (—180°) 


Determining Stability 
We will draw a Bode log-magnitude plot and then determine the value of gain that ensures that 
the magnitude is less than 0 dB (unity gai} at that frequency where the phase is 4 189°. 









A.J.A ۱ 52112012 Page 2 


Scanned with CamScanner 







valuating Gain and Phase Niargins 


Next we show how to eval : 9 
bg aluate the pain 

. 15 Irequenc . requency, Woy , wh th 1 
which is the - ت‎ G s at the magnitude plot to eme the "in — ۳ 1 
using the magnitude curve t re the magnitude curve to 0 dB. The phase margin is found by 
curve at that frequency. th o find the frequency, رون‎ Where the gain is 0 dB. On the phase 
y, the phase margin is the difference between the.phase value and 180*. , 


* Gain in: additi 
e "iter À— عدبم‎ gain that makes the system on the verge of instability. 
rgin: additional phase /ag that makes the system on the verge of instability. 


M (dB) 
Gain 
plot 
0 dB i | log @ 
1 
1 | CM 
1 
۱ 
Phase ۱ 
plot 
Phasc (degrecs} 
l 
يرح‎ ۱ 
180° log t» 


Oo, m cm 


Gain and phase margins on the Bode diagrams 


س صف wot‏ 


Je znow that iie charzcterislic equution zs 
| + KGH 0 
written for a particular gain. Then 
uS KGH = -1 

so that: a) when the amplitude of KGH is 1 or 0 dB and b) phase is 180° then we have the 
system just become unstable. Both the conditions have to be fulfilled. Hence, if we ask the 
question: How far is the system from instability ? We get two answers. 

@ The gain increment required at the frequency where phase is 180° and the phase increment 
required at the frequency where amplitude is 0 dB. The frequency where phase crosses 180° 
is called the phase crossover and the frequency where the amplitude crosses the 008 point is 





called the gain crossover. . 
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The figure below gives the bode piot for اک‎ 000 The point in the amplitude bode plot were 
gain crosses 008 is at frequency 6.08 rad/s. At this point the phase 1s -157.5 degrers In order 
to reach 180 degree I need a phase change of 22.5 degrees (-157.5 -( -180)) which 15 the 
Phase Margin. The phase becomes -180 degrees at 10 rad/s and here the amplitude is less 
than 1 or negative in dB. I have to raise the gain by 7.96 dB before the system becomes 
unstable. Hence, that is my Gain Margin. 


For the following transfer function find Gm and PM using MATLAB 


1000 


G(s) ˆ 5(s+5)(s+20) 


>> g= tf(1000, poly((O -5 -20[(( 


Transfer function: 
1000 


s^3 + 25 s^2 + 100 s 


۹ 
۱ 


له ٩۰‏ = .جب 2 - لسا ع وما 2 ده مة 


Bode Diagram 
Gm = 7.96 dB (at 10 rad/sec) , Pin= 22.5 deg (at 6.08 rad/sec) 


<< margin(g) e 











Meg ticle (E) 
$ 


g 
E 180—------ eed d 
-225 
-270 i E ۱ ۱ 
10 10° 10" 10° 10° 
Frequency (rad/s<c} 
>> [Gm,Pm,Ucg,Ucp] = margin(g) 
Gri = 2.5000 Pm = 22.5359 Weg = 10.0000 Wep = 6.0783 


>> Gm-207*10gi10 (Gm) 
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Feu the svete 


m shown in figure below, where Pee 


R(s) + جح‎ E(s) 3 1 tis) 
G(s) = CENE ہے‎ ۷ 1 G(s) | 
(s + 5)(s + 20) (s +50) | | 
8 


Do the following: 
a. Draw the Bode log-magnitude and phase plots. 


. b. Find the range of K for stability from your Bode plots. 


— ' e. Evaluate gain margin, phase margin, zero dB frequency, and 180? frequency from your Bode 
plots for K —210000. 


2, )وها‎ Me 


- —43-34 5 
13-99 3 ۱ 





Phase (degrees) 





Fras مخ‎ (rad/s) 


b. The phase angle is 180° at a frequency of 36.74 rad/s. At this frequency the 
gain is 99.67 dB. Therefore, 20log K = 99.67, or K = 96,270. We conclude that 


the system is stable for K « 96,270. 
c. For K = 10,000 , the magnitude plot is moved up 2010g10,000 = 80 dB. 


0 A 
Therefore, the gain margin is 99.67- 80 = 19:67 dB. The 180 frequency is 36.7 
3 | 
d/s. The gain curve crosses 0 dB at @= 7.74 rad/s, where the phase is 7 
rad/s. 


i 0 0 0 
We calculate the phase margin to be 180 —87.1 = 92.9. 
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MATLAB code 0 ظ‎ ۱ 
K= input ('Enter the value 512 ۰ 279 
G= tf(K,poly([-5 -20 -50])) $ G(s) H(s) 


rlocus (G) 

figure 

margin (G) —————— M" 
grid i | 

[Gn ,Pm ,Wcg,Wcpl-margin (G) 

Gm- 20*10gl0(Gm) ۶ gain margin db 


Enter the value of K :1 


K = 
x 


Transfer function: 


5*3 + 75 s^2 + 1350 s + 5000 


Gm = 9.6303e+004 
Pm = Inf 

Wcg = 36,752) 
Wcp = NaN 


Gm = 99.6728 db 
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K- 10000 


Transfer function: 


Gm = 9.6303 
Pm- 92.8979 
Weg- 36.7520 


- Wep= 7.7412 


Gm= 19.6728 db 











Magntude (dB) 





Q 
l ۳ تی‎ 
EX : Find the gain margin and phase margin if the value of K in each part the Problem ۰ ny 
5 V پ‎ 
3 S Q à 
K = 1000 , b. K = 100 and c. K - 0.1 Var) 
| 7 ۱ 
^ 
4, 
G Tat a 
“SO ل‎ C4 
-1 001 سي سهدت‎ CLE شی ای‎ ۶ af FH 
G.M. 53.6 dB + 
Freq: 6.63 rsd/sec F 
Stabia loop 
-15N 
0 =S 
= -90 TH F - # 1 ناه و ناه‎ P سح‎ 
E Hi ۱ 
& 0 سس را رہ ہے :ا‎ E MARLEY 
-270 سج 0111 ل اس ہہت‎ d Pii 
0 10 10° 10° | 10° 10° 
Frequency (rad/sec) 
K = 1000: 
‘For K = 1, phase response is 180 ° at w = 6.63 rad/s. Magnitude response is -53.6 dB at this 
frequency. ۱ 


For K = 1000, magnitude curve is raised by 60 dB yielding + 6.4 dB at 6.63 rad/s. Thus, the 
gain margin is - 6.4 dB. 


Phase margin: Raising the magnitude curve by 60 dB yields 0 dB at 9.07 rad/s, where the phase 
curve is 200.3" ۰ Hence, the phase margin is 180° - 200.3°= =- 20.3 ° 


K = 100: 


For K = 1, phas response i: 180? 2۶ o= £ 53 rad/s. Magnitude response is -52.5 dB at this frequency. 
For K = 100, magnitude curve is raised by 40 dB yielding — 13.6 dB at 6.63 rad/s. Thus, the gain 
margin is 13.6 dB. ۱ ۱ ۱ 

Phase margin: Raising the magnitude curve by 40 dB yields 0 dB at 2.54 rad/s, where the phase curve 
is 107.39. Hence, the phase margin is 1809-107.39 - 6 


K > 0.1: 


For K = 1, phase response is 180? at œ = 6.63 rad/s. Magnitude response is -53.6 dB at this frequency. 


For K = 0.1, magnitude curve is lowered by 20 dB yielding — 73.6 dB at 6.63 rad/s. Thus, the gain 
margin is 73.6 dB.. á 
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clase ail 
= input ('Enter the value of K d 


GH- tf(K,poly([-2 -4 -6])) $ G(s) H(s) 
rlocus (GH) 


figure 

margin (GH) 

grid 3 

[Gm ,Pm ,Wcg ,Wcp]=margin (GH) 

Gm= 20*1log10(Gm) ۶ gain margin db 


Enter the value of K :1 
K= 1 


Transfer function: 


1 


TL WT | |‏ ۱ 48 + و 44 + ۸2و 12 + 53و 


Gm = 480.0003 و‎ 


Pm= Inf 2 
Wcg= 3 P 
Wep= NaN 


Gm= 53.6248 db 
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For K= 1000 





Enter the value of K 0 ۱ 2 


K- 1000 


Transfer function: 


48 + و 44 + s^2‏ 12+ 53و 


Gm = 0.4800 

Pm = -20.2512 

Weg= 6.6333 

Wcp- 9.0654 
Gm = -6.3752 ۱ 
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Enter the value of K :0.1 
K= 0.1000 


Transfer function: 


553 + 12 s^2 + 44 و‎ + 48 


Gm = 0003 


Pm= Inf 
Weg = 6.6333 
Wep = NaN 


Gm = 73.6248 
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Use Bode plots to determine the rano; ee 
سه‎ ona the range of K within Which the unity feedback system 


K 
(s) (s+ 2)(s + 4)(s + 5) 


The low-frequency gain of G(s)H(s) is found b 
plot starts at K/40. For convenience, let K = 40 
each break frequency, 
the log- 


y setting 5 to zero. Thus, the Bode magnitude i3 
so that the log magnitude plot starts at 0 dB. At 7 
2,4, and 5, 820 dB/decade increase in negative slope is drawn, yielding ۱ 3 
magnitude plot as shown in figure below. The phase diagram begins at 0° until a x 
decade below the first break frequency of 2 rad/s. At 0.2 rad/s the curve decreases at a rate of 
-45° /decade, decreasing an additional 45°/decade at each subsequent frequency (0.4 and 0.5 


rad/s) a decade below each break. At a decade above each break frequency, the slopes are 
reduced by 45°/decade at each frequency. 


I 

/| LIN | 

$ 1 ۱ 

3 H ۱ 
۱ | 
T ۱ 

g ۱ 

= 





1 ۹ 
Frequency (rad/s) ۰ 


From the Bode log-magnitude plot must be less than unity when the Bode phase plot is 180°. 
Accordingly, we see that at a frequency of 7 rad/s, when the phase plot is -180°, the 
magnitude plot is -20 dB. Therefore, an increase in gain of +20 dB is possible before the — 
system becomes unstable. Since the gain plot was scaled for a gain of 40, +20 dB (a gain of 10) 
represents the required increase in gain above 40. Hence, the gain for instability is 

46 x 10 = 400. The final result is 0 < & > 400 for stability. 
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YQU 
NY QUIST STABILITY CRITERION: FUNDAMENTALS 
Thus far we have present ہے‎ d 
the Routh-Hurwits سا سیت کہ وب‎ 


the roots of the characteristic eq locus meth 


characteristic equation are al 
MATLAB. ilis 


g the stability of linear SISO systems: 
i od of determining stability by locating 
‘ation in the s-plane. Of course, if the coefficients of the 


Own, we can solve for the roots of the equation by use of 


Frequency-Domain Specifications 


1- Resonant Peak Mr The resonant peak 

Mr is the maximum value of] M(jw)}-In 

general, the magnitude of Mr gives ۷ از‎ 
indication on the relative stability of a 

stable closed-loop system. Normally, a M 
large Mr corresponds to a large maximum 
overshoot of thestep response. For most | 
control systems, it is generally accepted in 

practice that the desirable value of Mr 00 
should be between 1.1 and 1.5. 





2- The resonant frequency wr, 
is the frequency at which the peak 
resonance Mr occurs. 


=> 

ج 
us‏ 
= 
= 


3- Bandwidth ۷ ۱ 
The bandwidth BW is the frequency at | 
which|MGw)|drops to 70.7 % of, or 3 dB | 
down from, its zero-frequency value. In | 
general, the bandwidth of a control system 
gives indication on thc tı ansient-responsc 

properties in the time domain. A large 

bandwidth corresponds to a faster rise 

time, since higher-frequency signals are T 

more easily passed through the system. f 

Conversely, if the bandwidth is small, - (deg) 
only signals of relatively low frequencies 


are passed, and the time response will be Typical vain phase characteristics of a feedback control st 
slow and sluggish. Bandwidth also 
indicates the noise-filtering characteristics and the robustness of the system. The robustness 


represents a measure of the sensitivity of a system to parameter variations. A robust system 1S 
one that is insensitive to parameter variations. 
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۳ ines the stability.of a clo 
01110 :oraphical method that determines : 560 log 
, st criterion 1S a semigrap ۰ the 
i : b " vestigating the properties of the frequency -domain plot, fhe Nyquist plot, of the 
system by ‘in 2 5) ot Le) Specifically, the Nyquist plot o! L(S) is a plot of 


transfer function G(s)H( à 001 
LA "in the polar coordinates of Im[L(jw)] versus Re[L(jw)] as W varies from 0 to infinity, 


This is another example of using the properties of the loop transfer سر‎ to سب یسب‎ ۲ 
performance of the closed-loop system. The Nyquist criterion has the fol owing - res that 
make it an alternative method that is attractive for the analysis and design of control systems. 


1- In addition to providing the absolute stability, like the Routh-Hurwitz criterion, the Nyquist 
criterion also gives information on the relative stability of a stable system and the degr ee of 
instability of an unstable system. It also gives an indication of how the system stability may be 
improved, if needed. 


2. The Nyquist plot of G(s)H(s) or of L(s) is very easy to obtain, especially with the aid of a 
computer. 


3. The Nyquist plot of G(s)H(s) gives information on the frequency-domain characteristics 
such as Mr ,wr , BW, and others with ease. 


4. The Nyquist plot is useful for systems with pure time delay that cannot be treated with the 
Routh-Hurwitz criterion and are difficult to analyze with the root-locus method. 


EX 


Sketch the bode plot and Nyquist diagram for following system with unity feedback. 


cis) = 





1 
(s+2) 
S — jo 


1 
GU») = وص‎ 


1 2- jw 2- jo 2 W 
nh X بے سس‎ Ce = ca 
GQ ) 24jo 2-jo 440% 4407 4+ w? 











۱ ۱ ۱ ۱ ۱ 
Magnitude frequency response IG(jw)| = M(w) = I-A 

| | G(jo) = 0(w) = —tan"! =) 
Phase frequency response ۱ IG(jo) = 3 
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Nyquist Diagram 


Frequency (rad/sec) 


(Bep) ۵ 





(gp) ۵۸۸ 


Bode Diagram 


subplot (212) ;nyquist (G) 


subplot (211) 


bode (G) 


GÊ cf(1,[1 2]) 


~ 
r 


کہ 


-0.25 
Note: وا‎ = rad/s 


2 5 


31/3 


4 ris 
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RC 


Margin and Phase Margin via the Nyquist Diagram 


Im 











Unit circle 


ee ام سیت‎ gc 
e^ * 
5 * 


GH-plane 


Nyquist 
diagram 


Gain difference 
before instability 


Phase difference 
before instability 


Gain margin = G,,=20loga Phase margin = Dy = © 


Nyquist diagram showing gain and phuse margins 


EX 
Find the gain and phase margin for the following system if K=6, with unity feedback. | 


Ee ۱ أ‎ 
G(s)H(s) ` (5242542)(s+2) 





. ٦ 6 
G(jo)H (jw) = تی‎ + 25 + 2)6 + 2( 


نوزم و p‏ 0 6 ۱ 2 

a e CD j ہے‎ int) j £ 348 2 6[4(1 — a) — jo(6 — o*)] 

Vd S dd 16(1 — oy. + o?(6 — a? 
The Nyquist diagram crosses the real axis at a frequency of v6 rad/s. The real EP 
is calculated to be —0.3. Thus, the gain can be increased by 8ت‎ 3) = 3.33 before 
the real part becomes —1. Hence, the gain margin is 
Gy = 20 log 3.33 = 10.45 dB 

using Bode plots. G(jw)H(jo) has unity gain at a frequency of 1.253 rad/s 
At this frequency, the phase angle is —112.3° The difference between this 
angle and —180° is 67.7°, which is the phase margin. 
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Frequency (rad/sec) 
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Gain Margin (dB): 10.5 
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Bode Diagram 
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>> figure 
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EX Sketch the bode plot and Nyquist diagram for following system with unity feedback ۱۷ e / 


K =10,20,30,40,50,60,70,80,90,100 and check the stability in each case. 


K 
G(s) (s+4) 





for k= 10:10:100 

g- tf(k,[1 4]) 

bode +g) 

hold on 

end 

legend('k-10','k-20','k-30','k-40' ,' k-50' , 'k=60', 
k-70','k-80','k-90' , 'k-100') 1 
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fox k= 10 : 10 0 

g= ti(k,[1 4]) — 

nyquist (qg) a SD +1111 UTR 
hold on. 
end 

legend ('k=10' ',=99: 
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EX For the following system unity feedback. 
k 
G(s) ` (s+2)(s+4)(s+6) 


do the following: ۱ 
a. Plot the bode plot and Nyquist diagram. l u 
b. Use your Nyquist diagram to find the range of gain, K, for stability. 


c. Use your Nyquist diagram to find the gain margin and the | 80° frequency if k-100 . 
Bode Dlagram 
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b. Substituting jæ into G(s) = 


(48--1209?) — j(44@ — 03( ۱ 


simplifying, we obtain G(j@)= . The Nyqui 
piitying U®)= 5 S60" 17840 42304" Lhe Nyquist 


diagram crosses the real axis when the imaginary part of G(j@) is zero. Thus, the 
Nyquist diagram crosses the real axis at 0^ = 44, or @ = 444 = 6.63 rad/s. At 


Thus, the system is stable for K < 480.‏ سے 


this frequency G( jw) = 450 


If K = 100, the Nyquist diagram will intersect the real axis at 100/480. Thus, 


G&G, 20 log — 13.62 dB. the 180° frequency is 6.63 rad/s. 
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System k=1100 


Gain Margin (UB) -7.2 
At frequency (rad/sec): 6.63 
Closed Loop Stable? No 
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System: k=700 


` Gain Margin (dB): -3.28 
At frequency (rad/sec) 6.53 | 








Example of GM & PM 


٠ Second order system (5) = بع‎ 
Bode plot Nyquist plot 
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Comments 
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